COMBUSTION AND FLAME 97:393-417 (1994) 


593 


//•' - ' 

NASA-CR-202190 "V 

On the Influence of a Fuel Side Heat-Loss l “Soot”) Layer on a 

Planar Diffusion Flame 

INDREK S. WICHMAN* 

Laboratoire E.M2.C.. Ecole Central* de Pans. 92295 Chatenay-Malabry. France 


A model of the response of a diffusion flame (DF) to an adjacent heat loss or “soot" layer on the fuel side is 
investigated. The thermal influence of the "soot" or heat-loss layer on the DF occurs through the enthalpy 
sink it creates. A sink distribution in mixture-fraction space is employed to examine possible DF extinction, it 
is found that (i) the enthalpy sink (or "soot" layer) must touch the DF for radiation-induced quenching to 
occur, and (ii) for fuel-rich conditions extinction is possible only for a progressively narrower range of values 
of the characteristic heat-loss parameter. S R M R . Various interpretations of the model are discussed. .An 
attempt is made to place this work into the context created by previous experimental and computational 
studies. 


NOMENCLATURE 

a constant (see Eq. 13) describing reac- 
tion-zone heat loss (Eq. 14) 
a R same as a with radiation included 

A preexponential factor 

b reduced Damkbhler number, see Eq. 

14 

b e reduced Damkohler number at extinc- 
tion 

B E/3IRT,, see Eq. 26 

c an 0(1) constant in the definition of 77 , 

see Eq. 13. ii 

D, diffusion coefficients for species / 

D Damkohler number 

D r Damkohler number with radiant losses 

from the reaction zone 
DF “diffusion flame” 

e 2.718... 

E activation energy 

h enthalpy 

H “excess” enthalpy function, generally 

non-zero 

H a same as H with no losses; H 0 = 0 

HC "hydrocarbon” 

i bx blackbody spectral intensity given by 

Planck’s formula 
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k soot surface growth rate 

K function defined in Eq. 25 and evalu- 
ated in Sec. III.E 

L width of the physical domain 

L' density-coordinate reduced version of 

L , see Eq. 10 

Le, Lewis number of species i 

rh mass generated (destroyed) per unit 

volume per second (Sec. V) 

1 V “soot” number density, particles/ 

volume 

N r radiation number, see Eq. 2 et seq. 

P pressure 

P(x) polynomial defined by Eq. 21 
Q f heat released by combustion of fuel 

Q f nondimensional Q F , Q F = (1 + d>), see 

Eq. 2 et seq. 

Q r radiant heat flux 

Q r nondimensional radiant heat flux 

r ratio of reduced Damkohler numbers 

at extinction with and without radia- 
tion, r = b E R /b E . Also, reaction term 
in Eqs. 2 

31 universal gas constant 

S difference of temperature and enthalpy 

excess, S = r - H 
S? stagnation point 

time 

T temperature 

u velocity 

U a (x) heaviside step function U 0 (x ) 

_ / 0 . X < 0 \ 

\i. z > o j 
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creation/ destruction terms: phenomenology 

must therefore be abandoned and a set of 
kinetic equations for soot formation derived, 
along with a reasonable nucleation criterion. 

The approach used herein will be to exam- 
ine a physically simple model thoroughly. The 
emphasis will be placed on making mathemati- 
cally definite statements that can be translated 
into statements of physical fact, given the limi- 
tations of the model. We keep firmly in mind 
the principle that we are attempting to de- 
scribe — not simulate — the influence on DFs of 
radiant heat losses from particulates. 

We begin in Sec. II with the formulation of 
the problem including the heat loss function 
and the radiation term. In Sec. Ill we examine 
the results of the model predictions, including 
flame movement, extinction formulas, and 
bounds for radiant extinction. Then in Section 
IV we discuss the results, indicating points of 
strength and weakness in the model. The rela- 
tion of our work to previous studies is dis- 
cussed in Sec. V. Although literature reviews 
are usually placed in introductions we felt that 
for our model the literature review would be 
more useful if it placed the work in focus after 
the analysis was complete. Otherwise too many 
conceptual difficulties should confront the 
reader at the outset, serving only to obscure 
the subsequent deductions. We have at- 
tempted in most cases to examine the main 
features of the references we have cited. Fi- 
nally. Sec. VI presents a short set of point-form 
conclusions. 

II. FORMULATION 

I LA. Physical Discussion 

The following simplifications are employed in 
this study: geometrically, we consider the one- 
dimensional “stagnant film” diffusion flame 
(DF). The porous fuel wall, at temperature T 0 , 
L located at x — 0, parallel to the oxicizer waii' 

wtiil.il iiu. I * i g a..d Lb iJ^ulL.U til X “ U . 

The mass fractions of fuel and oxidizer at 
these walls are Y FF and Y 00 , respectively (see 
Fig. 1). Dynamically, we limit ourselves to the 
case of zero mean flow, so that the movement 
of species occurs strictly by diffusion. In addi- 
tion, we neglect the thermophoretic flow that 
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Diffusive riux or fuel Diffusive riux oi oxidizer 

Fig. 1. The physical configuration lor our model problem. 

is known to occur with particulates in regions 
of high thermal gradients, because we wish to 
examine the thermal and chemical influences 
of heat losses without factoring in dynamical 
complications that might make our efforts more 
difficult. In summary, we have a one-dimen- 
sional stagnant-film DF with no mean flow 
( u = 0) and no thermophoretic flow ( u T = 0). 
We also impose the steady-state condition 
{d(-)/dl = 0). For the combustion chemistry 
we assume that the reaction at the DF occurs 
through a single irreversible step, F + vO -» 
(l + u)P (on a mass basis), with high activa- 
tion energy. The "soot-formation mechanism,” 
which we hypothesize occurs on the fuel side 
of the DF, is assumed not to consume fuel. 
That is, only “trace" amounts of fuel are re- 
quired to make “soot” particulates. Also, since 
the “soot distribution” will be specified, we do 
not require a separate soot species equation. 
Nor is it necessary to consider a number-den- 
sity equation, since we assume that our “soot" 
particulates are simply a collection of immo- 
bile radiating masses located in a preassigned 
region on the fuel side of the DF. Strictly 
speaking, foei- :t " iv. '.ccd even f o discuss 
"pariicuiates” likcjso'?:” because none of the 
explicit features commonly associated with 
particulates appear in our analysis. The rele- 
vant features of our “soot particulate layer” 
are exclusively thermal', it produces only a re- 
gion of enthalpy loss that may alter the DF 
structure and cause extinction. 
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ti>) for the adiabatic flame temperature. We 
also observe that with heat losses the flame 
temperature will assuredly not rise to its theo- 
retical maximum even when — » x . Hence, we 

mav expect a slight redefinition later of the 
Damkohler number D in terms of a flame 
temperature somewhat lower than T,-. 

We note that the derivation of an excess-en- 
thalpy function may be achieved "physically. 
The conservation equation for the enthalpy 
takes the form of Eq. 4 when convective trans- 
port and body force effects and preferential 
species transport (Le, == 1) are_neglected. Since 
h = j/tfV; and h- t = h? + C p (T - 1°), it is 
easy to recover our nondimensional H. 

Finally, we observe that even in the most 
difficult and general case, such as when the 
radiation term depends on the spatial coordi- 
nate and the temperature and the fuel mass 
fraction, as long as no fuel — » soot depletion 
terms enter the species equations (Eqs. l.ii and 
I.iii) we can still define the mixture fraction 
variable Z = (4>y F 4 - 1 - yo)/(<£ + 1) that 
satisfies = 0 with Z = 0 at £ = 1 and Z = 

1 at £ = 0. This provides an important simpli- 
fication of the governing equations. The solu- 
tion for Z is Z = 1 — £, whereby 


(i) y 0 — (1 - Z) — (1 — Zf )( t H), 

(ii) » = Z - Z ; (r - H), 

( 5 ) 

where Zy — (1 + is the DF. location in 

the Z-coordinate system. Then the equation 
for r (the first of Eqs. 2) and the equation for 
H (Eq. 4) become 

(i) r zz = -(1 + < f>)Dr{H,T,Z ) ' 

+ N R (-dQ R /dZ), t(0) = t(1) = 0 ^ 
'iO %zz -<r {- dQ R /dZ). • ( 

H{0) = Hi 1) = 0. ' J 

( 6 ) 


Hence, the solution for r, y F , and y 0 is re- 
duced to the solution of two coupled nonlinear 
equations, Eqs. 6.i and 6.ii. Equations 5 and 6 


suggest defining 5 = r - H. giving 

ti) S zz = - 11 - o)Dr. 

5(0) = 5(1) = 0. I (7) 

lii) H zz = \ r R (-dQ R /dZ), 

H(0) = H(1)=0. 

where r(H . 5. Z) = (1 - Z - (1 - Z,)5] • 
[Z - Z f S]c\p[- (5i\ - 5 - H)/[ l - a(l - 5 
-H)}]. It is clear that some rather interesting 
behaviors may be expected, especially in the 
general case when the radiation term is a com- 
plicated function of Z. r, and perhaps other 
variables. However, we shall examine only the 
simple case when Q R is a prescribed function 
of Z. We shall see that even for this case many 
complexities arise. 


ILC. The Form of H(Z): 

The enthalpy defect H(Z) is obtained by inte- 
grating Eq. 4 or Eq. 7.ii twice. We consider the 
simple case when the radiant heat transfer 
term is a known, specified function of position. 
Then the integrations may be carried out ex- 
plicitly. Because of the eventual double inte- 
gration, we do not need to be particular in our 
choice for the radiant heat flux, dQ R (Z)/dZ. 
Hence, we let 


dQn 

dZ 


(Z) = U a (Z R .) - (/ 0 (Z R -), 


( 8 ) 


as shown in Fig. 2a. The quantities Z R - and 
Z R - are the boundaries of the heat loss zone. 
We note that Z R -> Z f and that there are no 
restrictions on Z R - other than Z R ~< 1, i.e., 
Z R - does not have to be “close” to Z R -. From 
Eq. 8 we see that dQ R /dZ is a “well” function, 
and that -dQ R /dZ is a “top-hat” profile. The 
solution for H{Z) (see Fig. 2b) is 


• | 0,- • 0 < Z < Z R -, 


i » ' 

* 


HfZ) - 


(Z Z R -){Z R ~- 


Z R -< z < z R ~, 
[0, z R ~< Z < 1. 



(9) 


Note that H z — 0 <*nd H — N R (Z R ~ 
Z R -) 2 / 8 at Z = (Z R -+ Z R -)/ 2, also shown in 
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This result, used in the definition of ,\' R alter 
Eqs. 2 and 3. gives 


,v» = 4| — 


PR 


Y, 


sR 


Y 


FF 


X 


*(V ~ T o) 


p 0 a 0 [Q F /(\ + <b)]/L* 


( 10 ) 


where L* — z 0 / p 0 is a length scale of the 
order of L, L*/L = L~ l f 0 L ( p/ p 0 ) dx, T f - T a 
= Q F Y FF /C p ( 1 + <b) and a 0 = \ 0 /p 0 C p . 
Hence. N R is the ratio of the blackbody radi- 
ant heat flux at T — T R to the flame heat flux, 
reduced by the length ratio L*/L RAD = 
L'k ? {T r , 1) < 0(1) and the ratio Y sR /Y FF . 

We shall see in the subsequent develop- 
ments that it is not N R by itself that is impor- 
tant. but the product of N R and the '‘soot 
layer” thickness A Z R in mixture-fraction space. 

Numerical estimates for N R can be made. 
We write p^/ p R = T R /T p and let T 0 = 300 K, 
T r = 2000 T K, where T is a nondimensional 
temperature. We also estimate Y sR /Y FF ~ 
1/20 and take Q F = 11355 cal/g (for 
methane). Then using p 0 = 3 x 10 3 g/cm 3 
and a 0 = 1.24 cm : /s we find N R = 
1.2L*(L*k p X\ + d>)T\ where L* is in cen- 
timeters. We expect L* ~ 2 cm and _L* K P ~ 
0(10" '), whereby N R — 0.14(1 + <£)T 5 . The 
factor 1 + 4> can range from 1 for d> < 0(1) to 
approximately 20 when_ is large, so that the 
quantity multiplying T 5 can range between 
0.25 and_approximately 5. A 10% increase of 
T r (i.e., f = 0.9) decreases N R by 40%, giving 
N r ~ 0.08(1 + & ), whereas a 20% decrease of 
T r decreases N R by about 70%. Consequently 
we expect N R to range from O(10" 3 ) to 
0(10" ‘) for small <£, to 0(1) - 0(10) values 
for large <b. 


ill. IlESULTS 

TIT. A. Temperature and Mass-Fraction Profiles 
in the Outer (Non-Flame) Regions 

On the oxidizer side y F = 0 in the lowest ap- 
proximation, giving y Q = \ — Z/Zf from Eqs. 
5.i and 5.ii and r = H + Z/Zf. Now if the 
“soot” layer is distinct from the reaction zone, 


i.e., if Z R -> Z.-. it has no influence on the 
temperature on the oxidizer side and r = Z/Z, 
=> dr/dZ = l /Z f there. In this case the oxi- 
dizer side is exactly the same as for a "non- 
sooting" flame. If. however, the leftmost edge 
of the soot layer and the flame zone slightly 
overlap, i.e.. Z R -< Z,-, then Z R -= Z f in Eq. 9 
(we anticipate that due to "oxidation” the dif- 
ference between Z R - and Z f will be small), 
giving 



'0. 0 < Z < Z R -= Z f , ^ 

i ~^{Z-Z f )(Z R .-Z),> 
[ Z R -<Z. 


and 




dr 

1 

dH 

1 iV R ±Z R 

dZ ~ 

Zf + 

dZ 

7 i 

Lf 


(11) 


as Z approaches Z f from the oxidizer side. 
Hence, the influence of the heat losses is to 
decrease the temperature gradient on the oxi- 
dizer side from its undisturbed value, 1/Z^. 
We observe that since d-r/dZ must be positive 
we obtain an upper bound for /V R AZ K , N R ±Z R 
< 2 /Z f = 2(1 + d>). This criterion may be in- 
terpreted as stating that when the heat losses 
become large the thickness AZ R can decrease: 
in other words, as N R increases we no longer 
need a thick “soot” layer to produce large heat 
losses from the flame. This argument can, of 
course, be generalized to parameters within 
N r . For instance, suppose that N R kZ R at ex- 
tinction has been determined and that A Z R is 
fixed, thereby fixing N R . Consequently if T R is 
changed at constant pressure, constant stoi- 
chiometric index d>, constant Q R , constant b , 
and so on, we must have x p (T R , l)Y sR T R = 
constant ; if k ? ~ constant, then a slight increase 
(or decrease ) of T R can be accompanied by a 
large decrease (or increase) of Y. R . Physically, as 
tne mean "souC i uyer tempti mure increases fewer 
soot particles are needed to produce the same 
heat loss. A 10% increase of T R , for example, 
decreases Y sR by 38%. 

We now examine the fuel side. Here we 
have y Q = 0 in the lowest approximation, giv- 
ing y F = (Z - Z / )/(l - Zf) and t = H + 1 
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toward the tuel side. We demonstrate this by 
writing the riame location as 

1 !\>'P dx 

l — (i - 0 1^ p dx 

This can be rearranged to yield 

,- x r ,-L 

I p dx = <b pax. 

f) J x, 

which we examine for the cases with and with- 
out heat losses. If the flame zone is negligibly 
thin we can write p = p 0 - ( p Q ~ P/XL - 
x)/{L - xf) on the oxidizer side, giving p(x f ) 
= p f and p(L) = p a and yielding 4>(L - 
x ■)( p 0 -r P f)/2 for the RHS of the above 
equation. For the LHS (the fuel side) we can 
use two different p distributions, the linear 
profile without heat losses,. p vl) = p 0 — (p 0 — 
p r Xx/xf, and a nonlinear profile with losses, 
p l2> = Po - ( p 0 - p f Xx/x f y-: p (2) is always 
larger than p^ l \ Substitution into the LHS 
above yields xf p a 4- pf)/- and xf p 0 + 2( p 0 
A pf)/ 2)/3, respectively, from which we find, 
after equating to the RHS, (.ty/L) (1) = d >/ (1 
- 1 - c b ) and ( x f /L) {1) = p4>/( P 0 /3 4- p(d> + 
2/3)), where p = ( p 0 + p f )P- We then form 
the difference. 


x F \ (1) l x f\ {2) ■ 

LI \ L ) ( 1 A <b ) ( 1 4- d> 4- M ) ' 



AJlid > o. 

Po + pf ) 


Hence, case 2 with heat losses ha's a smaller 
value of Xf/L, indicating that the flame is 
displaced to the fuel side. Although our 
demonstration has employed simple p distribu- 
tions, it is generally valid whenever p on the 
oxidizer side is unchanged by the heat losses 
and when the nondimensional reaction layer 
thickness, Ax^/L, is small. 

Because uur analysis is performed tn terms 
\jl iliC mixture tracnou z., iiic innuciivxb 
heat losses on Z f must be examined before 
making definite conclusions. In any case, the 
flame movement caused by density changes 
does not alter the flame location in the Z- 
coordinate, since Zf = (1 4- 4>)~ l is un- 
changed. 


III.C. The Chemical Reaction Zone 

First we analyze the case Z R -=Z ; -. so that 
H(Z.) = 0. We employ Eqs. 7.i, 9 with Z R - = 
Z f , 11. and 12. The use of the "S-equation" is 
of considerable benefit, because in all results 
except the presence of H in the exponential it 
resembles the standard DF “structure" equa- 
tion [1]. The gradients of S in the outer zones 
are identical to those for the zero heat-loss 
case, viz S z = (t — H) z = 1/Z^ on the fuel 
side (compare to Eq. 11 for t z ) and S z = (t — 
H) z = -1/(1 - Z f ) on the oxidizer side 
(compare to Eq. 12). We therefore define 
stretched variables. 


(i) 5 = 1 - (<J> + ar))/b(3, 

(ii) tj = J3(Z - Zf)c. 


as for standard DF asymptotic analysis [2]. We 
observe that S is presumed to approach unity 
in the reaction layer since 5 = r - H and H 
= 0 there, enabling 7 — » 1 as /3 — * °°. We shall 
subsequently see that when H(Z f ) # 0 we may 
still use Eq. 13.i because this implies r max = 1 
4- H{Zf) < 1, so that the flame temperature 
cannot attain its zero-heat-loss peak value. 
With the above substitutions the quantity in 
the argument of the exponential of Eq. 7,i 
becomes 


0(1 - S - H) 

1 - a(l - S - H) 


(-(O + a-n)/b + pH] 

a I d> 4- or) 

i-o —i P H 

P \ b 


which is in danger of being swamped by the 
heat loss term pH if appropriate measures are 
not taken. For Z near Zf, we write 

uh ( z, t ) 

H(Z) = rtlZf) - ( Z - t. ; d ~~ 

(Z-Zf) 1 d 2 H(Z f ) 

+ 2! dZ l + 

With H(Z f ) = 0, and dH(Z f )/dZ = 
NftAZf/2, and Eq. 13.ii, we find H(Z) = 
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When H(Z.) =f= 0. there is real overlap be- 
tween the flame and the ''soot" layer, so that 
Z R -< Z,. Nevertheless, the preceding analysis 
carries through in almost the same manner it 
we replace 0 in Eqs. 13. i and 13. ii by 0 R = 
/ 3T,/T, r . Here 7} R , the flame temperature with 
radiant heat losses, is smaller than 7}. It is 
obtained from the relation r R = 1 -4 H{Z , •), 
giving 7} R /7} = 1 + aH{Zf). The .Arrhenius 
exponent becomes 


0(1 - S - H) 

l - a(l - 5 - 77) 

= £ (1 ~ T *> 

1 - a(l - r R ) 


+ a R r > ] 

0 


-4 0( 0"‘). 

The term -0(1 - r R )/[ 1 - a(l - r R )] = 
-(£/lhXl/T /R - 1/7}) produces a modified 
Damkohler number, D R = D exp[ — (£/lX) • 
(1 /T fR - 1/7})], that is smaller than D. After 
some algebra we once again obtain the prob- 
lem given by Eqs. 14, with a R given by Eq. 15. 
a = 2Z f - 1, b/c = 2Z / (1 - Z f ), and b - 
[4(1 — <b)D R {Z f { 1 - Zy)/0 R ) 3 ] l/3 , in which D 
and 0 have been replaced by D R and 0 R . We 
note that since D R < D and 0 R > 0 the pres- 
ent reduced Damkohler number is always 
smaller than for the non-“sooting'’ flame. 


III.D. Extinction Criteria 

We deduce the criteria for the DF. with heat 
losses from the “soot” layer. There are at least 
two ways to do this. Both produce the same 
extinction zone on an N R hZ R vs. a (or Zf) 
plot, but the second method allows more physi- 
cal interpretation. 

In the first method we observe from Eq. 16 
that (i) - 1 < a R < + 1 and (ii) 0 < \a R \ < 1. 
Since a R = 2 Z, — 1 + (bf R &7-)Z f (l — Z ; ). 
we find tot criterion (i) after a short algebraic 
calculation that 1/2 < Z^ < 1 gives 0 < 
N r \Z r < 1/Zyr and 0 < Z f < 1 gives 0 < 
N r i\Z r < 2 /Z F . Thus, for the entire range of 
Z f , N R kZ R must be positive and smaller than 
2 / Z f . With criterion (ii) we are able to pro- 
duce a lower limit for iV R AZ R that exceeds 
when 0 < Z f < 1/2. We find the same 


limits for S R lZ R as in ii) when 1/2 < Z, < 1 
but when 0 < Z. < 1/2 we obtain 2(1 - 
2Z,)/(Z,(l - Z, )) < S R \ Z F < 2 /Z f . The 
extinction boundaries deduced here are shown 
in Fig. 4. They bound the lined region. The 
width of the extinction interval varies from 2 as 
Z,- — 0 to 4 at Zf = l /2 back to 2 as Z f -* 1. 
Observe also that extinction appears to be 
easier when 1/2 < Z, < 1 because the N R \Z R 
values required are fairly small. The physical 
reason for this is that the heat losses by con- 
duction from the reaction zone to the fuel side 
are greater than to the oxidizer side when 
1/2 < Z f < 1. We see this by writing Eq. 16 as 
1 - a = 2/(1 4- y), where y = |r z | fue i/|r z | 0 , ; 
hence, y > 1 when 0 < a < 1. When the "soot” 
layer is on the fuel side the heat losses to it are 
amplified, making extinction easier. When 0 < 
Z f < 1/2 these two factors act in opposition, 
making extinction more difficult. 

In the second method we examine the ex- 
tinction criteria for Eq. 14 in detail. It is well 
known that for a fixed value of a R , Eq. 14 
behaves as follows: when b > b £ , where b E is 
the ‘"extinction” Damkohler number, there are 
two solutions, only one of which we can accept 
as “physically realistic”. When b < b E there 
are no solutions. Each value of a produces a 
single value of 6, enabling the determination 



Fig. 4. jVfl AZfl versus <x and When 7 ^ 7^ E 
extinction can occur only in the lined region. The dashed 
vertical lines at Z f - 0.08 and Z f - 0.4 indicate bounds 
proposed by Kennedy et al. [14], outside of which they 
assert that the soot growth rate is zero. When Tf*.E>T,.E 
extinction must occur in the cross-hatched region. 
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At extinction we can equate our two ratios 
r- to obtain 


/ ( 1 t! Ri) 

/(I — m) 


T, 


fR. £ 


-E 


T f.E 


exp 


3917V r 


T, 


. E 


7V 


- 1 


f. £ \ 1 fR. £ 


(. 22 ) 


There are two cases, T fR E < T f E and T fR E > 
T fE . 

Case (i) T fRE < T f E 

Here we have r E < 1. This can be satisfied only 
when 1 — |a R | < 1 - 14, or |a R l > !4 (see Fig. 
5). We recall from Eq. 15 that a R > a is always 
true. Hence, extinction can occur only when 
the two inequalities a R > a and In,/ > 14 are 
satisfied. When a > 0 these inequalities are 
obviously satisfied and produce the admissible 
range of values 0 < N R EZ R < 2/Zf = 4/(1 + 
a). When a < 0 (or 0 < Z f < 1/2) we always 
have a R > a, but \a R \ > 14 requires \a R \ - a -f 
(,V ft AZ fi )Z,(l - Z f ) > 14 = -a, whereby 

N r ±Z r > 8(-a)/(l~a-) = 2(1-2 Z f )/ 

{Z f {\ - Z f )), which is the same set of limits 
deduced before. In summary, 


Case tii): T, R E > T, E 

Here we have r r > l. so that a R > l and 

< 14. These conditions restrict us to - l < u 

< 0 (0 < Z, < 1/2). where we find the lower 
limit ,\/A Z R = 0 and the upper limit \ R EZ H 
= 2(1 - 2Z,)/[Z.(l - Z. )l, which coincides 
with the lower limit in Eq. 23 for the case 
— I < a < 0. This extinction region is shown in 
Fig. 4 as the cross-hatched zone. 

We therefore have demonstrated that this 
approach allows us to describe the region be- 
low N r EZ r = 2 /Zf. 

Our results raise two questions. (1) Does the 
region produced by Tf R E > Tf E make 'physi- 
cal sense"? This is debatable, because for con- 
ditions other than extinction Eq. 20 demon- 
strates that r(Z f ) = b R /b is always less than 
unity. How then is it possible for b R /b at 
extinction to suddenly rise to above-unity val- 
ues? We suspect that this region is indeed 
nonphysical. (2) Can fuel-poor DFs with 6 < 1 
(or 1/2 < Zf < l, or 0 < a < 1) even produce 
enough ‘Noot" to self-extinguish? According to 
our model they can indeed. In practice, how- 
ever, the condition 6 < 1 places tight restric- 
tions on possible extinction conditions. For the 
general hydrocarbon-oxygen reaction 


f Z(1 - 2 Z f ) 8(-a) 


Z,(l - Z f ) l-fl- 


jV r AZ r > \ 


0 <z,< 


0 ; 




and 


N r \Z r < — 


(23) 


n v z, < i (-i 


<-iv 


(24) 


This produces the same lined region of Fig. 4 
as previously, but with the added restriction 

TfR £ < T f. £ • Here the DF with heat losses 
extinguishes at a lower flame temperature than 

the DF without heat losses. 


C„H m + {n - r m/ 4)(0 : - -tN ; ) 
-» nCO z -r (m/2)H : 0 
4- (n -r m/4)xN : , 


we have v— 32(n 4- m/4)/(12n + rn). which 
ranges from 3.08 to 4.23 when m/n ranges 
from 1 to 5. Hence, d> < 1 implies T E{: /Yqq < 
\/v, which ranges from 1/4 to 1/3. For air 
(Y 00 = 0.23) this means Y FF e (0.055.0.075). 
which is a rather narrow interval, possibly be- 
low the extinction limit and certainly so fuel 
poor that expectations of “<oot” production 
may be in vain. For' eievated-oxygen environ- 
ments perhaps encugn boot may torm to 
weaken the DF. 

On the fuel-rich side of Fig. 5 the •'physi- 
cally realistic” interval for extinction (the lined 
region in Fig. 4) becomes progressively more 
difficult to attain as Z ; decreases. For HC-air 
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heavier than in (b). 


0.2 to 1.0 in increments of 0.2. Thus, Z^ ranges 
from 0.2 to 0.55. The abscissa is A R AZ R - 
0V R AZ R ) min , which has a 'range of 2.5 when 
Z, = 0.2, a range of 4.0 when = 0.5 and a 
range of only 3.6 when Zj — 0.55; this is why 
the trend to increased range of the successive 
curves for Y QO = 0.2, 0.4, 0.6, and 0.8 is broken 
when Y 00 = 1.0. In Fig. 6b we have Y FF = 0.6 
with Y 00 = (0.2, 0.4, 0.6, 0.8, 1.0). Here the 
range of N r i\Z r continually increases. The 
same behavior is found in Fig. 6c, where Y FF — 
1.0. The data of Figs. 6a-6c are normalized by 
dividing by the maximum value of each final 
data point for each curve. This normalization 
improves the correlation of the data as Y FF 
increases, see Fig. 6d. More sophisticated nor- 
malizations are clearly possible, but this is an 
exercise we do not undertake here. 


TABLE i 


Parameter Values Used to Generate Figs. 6 and 7 


Parameter 

Numerical Value 

E 

29.100 cal/mol 

c. 

0.323 cal/g-K 

Qf 

1 1.355 cal/g (fuel) 

V 

-i.O 

To 

2S9 K 


It is perhaps more useful to focus on A = 
a\H{Z f )\ since this is the temperature decre- 
ment produced by the heat-loss zone. For the 
same Y FF and parameter values, the results for 
A vs. jV r AZ r - (N R ^Z R ) mln are shown in Figs. 
7a-7c. Observe the similarity to Figs. 6a-6c. 
The correlation analogous to Fig. 6d is shown 
in Fig. 7d. 

IV. DISCUSSION OF RESULTS 

We have made numerous simplifications dur- 
ing our analysis. It is necessary to examine 
these in order that future studies might intro- 
duce the necessary refinements. 

The radiation term in Eq. 1 (or Eq. 4, or Eq. 
7.ii) was simplified by ignoring its functional 
dependence on r, y F , etc. It was replaced by 
the simple model function of Fig. 2b that al- 
lows decoupling of Eq. 7.i and 7.ii and facili- 
tated the subsequent analysis. Later, in Sec. 
II.D, we assumed the soot layer was an opti- 
cally thin medium, which allowed us to relate 
dQ R /dZ in Eq. 7.ii to the soot mass fraction y s 




Fig. 7. (a) A = a\H(Zf)\ versus N*AZ S - (N*AZ R ) raiB . 
Same conditions as Fig. 6a. 
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limit when k, 1 d(-)/dx <k: I. The advantage 
of this formulation is that the integrodifferen- 
tiai equation of radiant transfer is replaced by 
a nominally tractable ordinary differential 
equation. This might enable analytical solu- 
tions to be found, or. at least, numerical results 
to be obtained without employing sophisticated 
and expensive computational software and 
hardware. Just as importantly, it preserves the 
modeling imperative, the notion that in model- 
ing we wish to describe what happens, not 
simulate it. 

From these discussions of the purely radia- 
tive features we see that much "structure for 
S and H remains to be uncovered from Eqs. 

7.i and 7.ii. We have, through the specification 
of H(Z), produced a simplified theory that 
only touches the surface of this ‘‘structure. 

The approximation made for the reaction 
chemistry was standard, in contrast with stud- 
ies that include detailed chemistry. In order to 
produce a model for the chemical reaction 
processes that describes fuel breakdown and 
soot formation, however, we do not require a 
detailed numerical simulation, unless we are 
certain that soot production is initiated upon 
the appearance of a specific "precursor. All 
we require is a qualitative reaction sequence 
that resembles the actual sequence in its most 
important features. With our fast-chemistry 
model, we determined that the "soot layer 
had to touch the reaction zone to influence it. 
This conclusion may change when more than 
trace amounts of fuel are required to produce 
“soot”, but the implications are clear: the de- 
pendence in the reaction zone on T through 
the Anhenius factor exp(-£/9UT) is strong 
enough to drive the chemical reaction to com- 
pletion even when broad regions of heat loss 
exist nearby. Only when the losses severely 
interrupt the zone of maximum reaction rate, 
which for <t> < 1 (fuel-poor or oxvgen-rich) oc- 
curs slightly to the oxidizer side of Z f , and for 
<i> > U occurs slightly to tne fuel side of Zf, 
does tne possioitity of uF extinction exist.- 
The extinction mechanism is thermal be- 
cause the sole influence of “soot” is to create a 
heat-loss zone. In fact, there is no reason to 
presume that our virtual radiating particulates 
have their origins in soot. The heat-loss zone is 
generic and may perhaps be best thought of as 


follows: imagine a collection of small chemi- 
cally inert spherules distributed on the fuel 
side of the domain, because oxidation destroys 
them when they penetrate the flame. Above a 
certain temperature they radiate energy out of 
the domain 0 < .v < L. We may assume that 
the heat flux leaving each point of the heat-loss 
zone is identical, giving for - dQ R /dZ the 
"top-hat” profile already discussed in Sec. II. D. 
This amounts to defining a characteristic loss- 
zone temperature T = T R , viz, T R (T R — T 0 ) 

= - T*)dZ\ the constancy of 

Y s on the fuel side would eliminate the need 
for a "soot” mass profile like the one described 
by the “top hat” profile of Fig. 2a. The result 
would be a radiant heat-loss zone between two 
temperatures (the minimum, T Smin and the 
maximum. 7) R ) separated by the distance AZ R 
in mixture fraction space. 

In this article we have moved the "soot 
layer” at will or, reverting to the above intepre- 
tation, we have defined the bounding tempera- 
tures of the heat-loss layer at will. Heat-loss 
regions in real DFs, however, usually locate 
themselves through a web of interactions in- 
volving nucleation. particle growth, convective 
flow, thermophozetic flow, etc. If we presently 
ignore the details of nucleation and panicle 
growth and assume that mature panicles spon- 
taneously appear, we must still account for the 
imbalance of physical forces. In our study we 
have neglected convection and thermophore- 
sis. It is known that convective flows of fuel 
toward the DF and thermophoretically induced 
paniculate motion away from the DF strike a 
balance that dictates the soot layer location. 
Oddly enough, without convection but with 
thermophoresis, the stable condition for our 
problem is a steady thermophoretic flow of 
particulates toward the fuel wall. Thus, the 
entire fuel side would be populated with car- 
bon particulates streaming steadily towards the 
fuel wall. This suggests locating the soot zone 
according to strictly thermal crnenc. which in 
fact appears ic de tne most p/iyswally jetj-con- 
sistent way to interpret the model examined in this 
article. 

For the flow field there are many complica- 
tions. In this study they have been eliminated 
by ignoring the flow, though it would be rela- 
tively straightforward to introduce a convective 
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portant: increase it to form more soot, de- 
crease it to produce less soot. 

Kent et al. [9], with their Wolfhard- Parker 
burner, found that the soot particulates are 
generated very close to the reaction zone, that 
the soot formation rate has its peak values 
about 2-3 mm from the reaction zone on the 
fuel side, that the soot volume fraction and 
particle diameter profiles peak at about 5 mm 
on the fuel side of the 7-peak, that near the 
flame base the soot formation rate becomes 
negative (oxidation), and that the soot particu- 
lates can be convected through the DF. 

The 1985 study by Smyth et al. [10] has an 
extensive bibliography (59 references) and cer- 
tain succinct and direct conclusions, one of 
which we now quote. With their Wolfhard- 
Parker burner for CH 4 /air DFs these authors 
determine that '‘...condensation reactions are 
found in a localized region approximately 2-3 
mm on the fuel side from the zone of maxi- 
mum temperature. These intermediate FICs 
exhibit maximum concentrations at tempera- 
tures of approximately 1300-1650 K. The earli- 
est soot particles are detected at the high 
temperature edge of this region, which is rich 
in aromatic compounds and other unsaturated 
HCs.” The Kent et al. [9] and Smyth et al. [10] 
studies are in general agreement over the prox- 
imity of the soot layer to the reaction zone. 
But Wolfhard-Parker burners have substantial 
rates for free and forced convection, perhaps 
negating the thermophoretic flows and produc- 
ing a bias towards proximate soot layers and 
DFs. 

The studies of Santoro and Colleagues ex- 
amined soot growth along individual particle 
paths. They used a coannular DF and per- 
formed extensive particle, temperature and 
velocity-field measurements using laser ve- 
locimetrv (LV) [11], They estimate a maximum 
thermophoretic velocity u T max of 4 cm/s and 
observe that in the region where u T — u Tmix 
the radial velocity ranges Frrr.cen iu ana .20 
cm/: forced and free convection dom- 

inate thermophoresis. Soot is formed in the 
annular region of the flame, where the soot 
volume fractions are maximum, 3> mM . As the 
burner flow rate increases the path line for <t> 
“shifts” outward to larger radial values. The 
velocity field is dominated by buoyancy, mak- 


ing it insensitive to changes in the fuel flow 
rate, but the soot field is strongly dependent on 
t, through the fuel flow rate. The two major 
contributions to increased soot formation are 
increased t r and changes in the flame shape 
caused by the ''shift” of Increases in the 
fuel flow rate do not increase soot formation 
rates. The values of are 0(10) to 0(100) 
times larger than near the centerline of the 
flame axis, except for the region close to the 
flame tip, suggesting that soot is indeed able to 
penetrate (diffuse?) to the center of the flame 
cone given sufficient time. Unfortunately, there 
is no information on the figures of the relative 
locations of the soot, flame, and oxidation re- 
gions. Nevertheless, it is clear that increased 
fuel flow rates enable soot particles to move on 
path lines with larger radial values, indicating 
that regions of soot formation respond to a 
change in flow rates. 

A recent study that supports this statement 
is that of Jackson et al. [12], for /xg droplet 
burning of heptane and heptane/monochloro- 
alkane mixtures. The latter fuel soots more 
than the former because the Cl + radical 
“scavenges” the OH“ atoms that usually in- 
hibit soot formation. These authors suggest 
that the spherical “soot shell” is relatively close 
to the droplet surface. Hence, contrary to the 
previous studies [8-11] it appears the inwardly 
directed thermophoretic forces can compete 
with the outward convective flow. They indi- 
cate that “soot shell” formation occurs where 
F d + F t - 0, where particle drag and ther- 
mophretic forces cancel, that the force balance 
depends on agglomerate size and gas proper- 
ties (mostly the viscosity p. and the diffusion 
coefficient D), and that the soot formation can 
be minimized by ensuring the formation of 
smaller particulates that do not break through 
the DF, but are consumed by it. This study 
raises questions about the thermophoresis- 
convection competition. The force balance sug- 
gests tliai ihermopnoresis can be very impor- 
tant. if soot pztticuiaic.i nucleate at icaipcta- 
tures above 1650 K (as suggested in [10]), there 
must be a rather vigorous migration from the 
high-T near-DF region to the low-7 near- 
surface region. Can such rapid migration away 
from the regions abundant in aromatic com- 
pounds and unsaturated HCs [10] still allow 
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not include 7 in their calculation of d 5 . so that 
although 7 drops once the flame is crossed 
from the fuel to the oxidizer side, the corre- 
sponding drop in is much smaller than it 
might be. .An ad hoc attempt to include the 7 
effect (which they say simulates the tempera- 
ture drop due to radiation) improves the 
agreement only slightly. They also find that d> 
depends strongly on k. increasing by a factor 
of 10 when k increases only by a factor of 2.5. 
They also note that k is actually a function of 
Z and time r. and that k decreases as t in- 
creases. because "young'’ soot has a higher 
surface growth rate than "old” soot. They con- 
clude by drawing attention to the need for 
accounting more precisely for the influences of 
radiation on the 7 distribution in sooting 
flames. A means must be found to introduce 
variable 7 into the equations for d>, N, etc., 
and radiative loss terms as functions of 
d>. N , ... must also be derived. 

The study of Pagni and Okoh [15] examines 
certain fundamental questions of soot-flame 
interaction using a two-dimensional laminar 
boundary-layer flame on a pyrolyzing slab. This 
is really a multidimensional version of Fig. 1. 
Perhaps because of the strong convective flows 
they find negligible thermophoretic influences. 
They employ a simplified phenomenological 
expression of Arrhenius form to describe the 
soot growth rate. Based on existing empirical 
data they state that the "... optically thin ap- 
proximation is valid for incorporating panicle 
radiation in the gas boundary layer equations.” 
The previous study with perhaps the most 
relevance to our work is that of Leung et al. 
[16]. It advances a mechanism for soot forma- 
tion in DFs that with refinement should be 
able to eventually absorb the phenomenologi- 
cal models. The authors begin by making sev- 
eral observations: (i) “ . . . measurements indi- 
cate that soot formation is dependent on the 
breakdown path of the fuel and the presence 
of pvrolysis products Mien as acetylene., and 
pclyans«iuraieu cvciical HO such as 
benzene ...” (ii) pyrolysis products "tend to 
show similar profiles with different 
magnitudes . . . so that numerous “critical 
species” exist, like C 2 H 2 , C 6 H 6 , C 4 H 6 , etc.; 
and (iii) the detailed models are too complex 
to use in “flows of practical importance for 


vears to come” and the global (phenomeno- 
logical) models are of dubious validity "under 
conditions different from those under which 
they were originally formulated.” They con- 
struct a soot formation mechanism for DFs 
containing five components. The first compo- 
nent is the chemical reaction mechanism that 
predicts at least qualitatively correct levels of 
C-.H ; , the "precursor” they define as the con- 
duit to sooting. The second component is the 
formation of incipient particles in regions of 
high precursor concentration. Although some- 
thing less than 10% of the total soot mass 
consists of these incipient particles, they are 
important because they provide the initial con- 
ditions for the subsequent growth phase. The 
second step, the “nucleation rate,” is assumed 
proportional to the precursor concentration. 
C c , H .. The minimum panicle size is taken to 
be 100 C atoms. The third step is surface 
growth by adsorption. For their model they 
employ the results of previous experimental 
studies [17-19]. The fourth step is soot oxida- 
tion to CO, which the authors assert occurs in 
a narrow region close to the flame on the 
oxidizer side. They also state that their 
“ . . . predictions are comparatively insensitive 
to this reaction step.” They employ 0 ; as the 
oxidizer, though strong arguments have been 
made for OH. The fifth component of their 
model is a particle agglomeration step. Inter- 
estingly, they employ detailed kinetic equations 
that have been derived for adiabatic flames: 
hence, their detailed mechanism should be ap- 
plicable for lightly sooting flames. For the 
model equations only thermophoretic trans- 
port is specified for the soot, and a rather 
simple ad hoc correction was used to simulate 
the radiant heat loss from the flame: for the 
temperature distribution they wrote 7 = r ad d 
- C(7 ad /7 max ) 4 ), where 7 ad is the local tem- 
perature calculated adiabatically and 7 max is 
the maximum value of 7 ad . Obviously this cor- 
rrr'ijn cannot change the qualitative shape of the 
;?mp-~-;ure curve, it can* tOi add an inflection 
point in regions of high heat loss, for instance, as 
we have in Figs. 3a and 3b. Nevertheless, d> is 
fairly well predicted, with discrepancies at- 
tributed to “uncertainties in the gas-phase 
chemistry and the simplified treatment of the 
nonadiabaticity ...” They conclude that for the 
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strong temperature gradients in the particu- 
late layer. 

7. Although our model is a highly simplified 
version of any available physical (experi- 
mental) counterpart, our discussion of the 
literature shows that many of our simplifi- 
cations. restrictions and assumptions are 
consistent with empirical observations. Plac- 
ing the “soot” layer adjacent to the reaction 
zone is in reasonable accord with the evi- 
dence. It is likely that (convection) » 
(thermophoresis) in the studies we have re- 
viewed, except for Ref. 12. 

8. Arguments were presented supporting the 
need for a reliable and physically reason- 
able physicochemical mechanism for soot 
formation. The phenomenological ap- 
proach, with its strong intuitive appeal, re- 
quires the existence of a convective flow. 
Although our zero-convective limit may not 
be realistic, it focuses attention on the defi- 
ciencies of phenomenological modeling. 

9. An interesting set of equations (Eqs. 7.i, 

7.ii) was obtained for describing the interac- 
tion of a DF and an enthalpy-loss region. 
When H = only these two equa- 

tions need be examined. When the radiation 
term dQ R /dZ also depends on the soot 
mass fraction, which itself depends at the 
very least on the fuel mass fraction and the 
temperature, we have a coupled system of 
three nonlinear ordinary differential equa- 
tions. In the general case the situation is 
even more complicated, for the soot mass 
fraction depends on the “initial conditions” 
provided by the pyrolysis chemistry, the nu- 
cleation sequence, the surface growth rate, 
etc. As always, the number of coupled non- 
linear equations increases. 

The author gratefully acknowledges a stimulat- 
ing discussion in the winter of 1991 with Dr. 
.Michael Deiichatsics, ns well as me receipt of an 
irueicsti/ig unpublished rzputt fran hu n._ i nc.se 
together suggested that theoretical models could 
indeed be developed to profitably examine this 
complex problem. This research was supported by 
NASA under the micrograuity combustion pro- 
gram, contract number NAG3-1271, monitored 
by Dr. K. Sacksteder. 


REFERENCES 

1. Linan. A.. Acta Astronaut.. /: 1007-1039 ( I9~4l. 

2. Williams. F. A. Comoustion Theory, 2nd cd.. Ben- 
jamin-Cummings. CA l 1 935). 

3. Tzeng, L*S.. Atreya. A., and Wichman. I. S.. Com- 
bust. Flame SO :9a ( 1990). 

A Daguse, T. Sounam. A.. Darabiha. N., and Rolon. 

J. C.. IV Internationai Seminar on Flame Structure, 
18-22 August, Sovosibtrsk. Russia . 1992. 

5. Siegel. R.. and Howell. J. R.. Thermal Radiation Heat 
Transfer . McGraw-Hill. New York, 1972. 

6. Osizik. M. N.. Radiative Transfer and Interactions with 
Conduction and Convection, Wiley, NY, 1972. 

7. Joulin, G.. and Deshaies. B.. Combust. Sci. Technol. 
47: 299-315 (1986). 

8. Wagner, H. Gg., Seventeenth Symposium (Interna- 
tional) on Combustion. The Combustion Institute, 
Pittsburgh, (1979). pp. 3-19. 

9. Kent, J. H., Jander. H.. and Wagner, H. Gg., Eigh- 
teenth Symposium (International) on Combustion , The 
Combustion Institute. Pittsburgh, (1981). pp. 
1117-1126. 

10. Symth. K. C., Miller. H., Dorfman. Mallard. W. G. 
and Santoro, R. H., Combust. Flame 62:157-131 
(1985). 

11. Santoro, R. J., Yeh, T. T., Horvath, J. J., and Semer* 
jian, H. G., Combust. Sci. Technol. 53:89-115 (1987). 

12. Jackson. G. S., Avedisian. C. T., and Yang, J. C., Int. 
J . Heat Mass Trans. 35:2017-2033 (1992). 

13. Kent, J. H.. and Honnery, D. R., Combust. Flame 
79:287-298 (1990). 

14. Kennedy, I. M., Koiiman. W., and Chen, J.*Y.. Com- 
bust. Flame 81:13-85 (1990). 

15. Pagni, P. J., and Okoh. C. I.. Twentieth Symposium 
(International) on Combustion , The Combustion Insti- 
tute. Pittsburgh, 1984. pp. 1045-1054. 

16. Leung, K- M., Lindstedt, R. P., and Jones. W. P., 
Combust. Flame 87: 289-305 (1991). 

17. Harris, S. J., and Wiener. A. M., Combust. Sci. Tech - 
nol. J7 : 155— 167 (1983). 

18. Harris. S. J., and Wiener. A. M., Combust . Sci. Tech - 
nol. 32:267-275 (1983). 

19. Vandsburger, U., Kennedy, I. M.. and Glassman. L 
Combust. Sci . Technol. 39: 263-285 (1984). 

20. Sohrab, S. H.. Linan. A„ and Williams, F. A., Com- 
bust. Sci . Technol. 23:143-154 (1982). 

21. Bonne, U., Combust. Flame / 6 : 147— 159 (1971). 

22. Hottel, H.. and Saronm. A. F., Radiation Transfer , 
McGraw-Hill. New York, 1967. 


APPENDIX 

Here we generalize the formulation of Sec. II 
vis a vis H{Z). First, however, we reexamine 
the functional form we used there in slightly 
greater detail. We observe the following: (i) 
The derivative of H(Z) is discontinuous at 
Z R - and Z R ... with values ±N R SZ R /2 inside 
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If K/V', does not vanish. H zz is positive 
everywhere on the interval 0 < Z < 1. It is zero 
only at the endpoints. Hence. H topologically 
resembles Fig. 2b with Z R - = 0 and Z R - = l. In 
reality. Y,/Y xR will vanish at least on the side 
Z < Z. due to oxidation. If it vanishes abruptly, 
as in Fig. 2b. no further discussion is necessary': 
but if it vanishes smoothly, producing a "dos 
ear” near Z = Z,- we may argue that the oxi- 
dization of the soot produces a local heat 


source. On the fuel side we do not expect anv 
such heat sources: hence, either Y, vanishes 
abruptly or it remains non-zero to the wall, 
where the condition 1=7, puts an end to the 
region H zz > 0. In either case. H zz should 
approach zero with a nonzero slope. This is 
qualitatively identical to what we have already 
done in Fig. 2. 
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